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Search for tt t^Vt with hadronic tagging at Belle 
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We search for the process B^ —>■ using the full Belle data set of 711 fb“^, corresponding 

to 772 X 10®-SB pairs, collected at the T(4S) resonance with the Belle detector at the KEKB 
asymmetric-energy e''"e“ collider. We reconstruct one B meson in a hadronic decay and search for 
the B® —>■ process in the remainder of the event. No significant signal is observed and an 

upper limit of B (B® —>■ 7r“r"'"i/r) < 2.5 x 10“"* is obtained at the 90% confidence level. 

PACS numbers: 13.20.He,12.60.Fr 


I. INTRODUCTION 


The decay B® —>■ [T] is mediated by the W'^ 

boson via the b ^ u transition. The transition amplitude 
is described by [2] 


(tt = /+(g^) 






On 




-9m. 


( 1 ) 


The form factors f'^ and /° can be computed from 
QCD light-cone sum rules [H |3] for < 16 GeV^/c^ 
and lattice QCD [3HS] for > 16 GeV^/c'^. Various 
parametrizations exist to interpolate between the two 
regions. In this study, we use the parametrization in¬ 
troduced by Bourrely, Caprini, and Lellouch (BCL) [7], 
which can describe both form factors in < 

{rriB — itLtt) . The parameter values are taken from 
Ref. [2]. 

It has been stated Hi that the differential ratio 

dr (B — TTTVr) 
dr (B —)• 7r£z/^) /dq^ ’ 


with p and q being the momentum transfers to the pion 
and lepton pair, respectively. 


i = e, fj, 


( 2 ) 
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can be used as a test for the Standard Model (SM) as it 
depends solely on the ratio of the scalar and vector form 
factors ■ The CKM matrix [3] element \Vub\ enters 

both differential branching fractions and cancels in the 
ratio. 

In new physics models like the two-Higgs-doublet 
model (2HDM) (TUI fTTj . the decay —)■ tt can 

also be mediated by a charged Higgs boson. Possible 
contribution of a H'^ and other couplings in the 2HDM 
and MSSM [T^ [T3], which would affect the branching 
fraction and the differential ratio of branching fractions, 
have been evaluated in Refs. [2] and [HHiT]. 

The decay —)■ 'K~T'^Vr has not been ob¬ 

served, nor has an upper limit on the branching 
fraction been obtained. Recent results [5] on the 
two form factors obtained from a joint fit to (2-1-1)- 
flavor lattice QCD calculations and B —>■ Triv data 
from Belle m US] and BaBar na nn result in 
B (5° —)■ 7r“T+i/7-) jB (R° —)■ 7r“£+i/f) = 0.641(17) and 
B (B° ^ TT-T+vJ) = 9.35(38) x IQ-’^ [H]. 

The signal decay is reconstructed in the four one-prong 
decays of the r lepton, t~ —>■ with £ = e or 

T~ and t~ P~Vt, corresponding to 72% 

of all r decays |13]. The most powerful decay modes 
are the two aforementioned hadronic r decays and the 
T~ — > e~VeMr mode, while the t~ —> decay mode 

does not improve the final expected significance. This is 
mainly due to low muon momenta in the signal decay 
and the resulting low muon identihcation efficiency. The 
result of this analysis is based on the three most powerful 
T decay modes. 


II. DATA SAMPLE 

The search for B^ —>■ described in this paper is 

performed on the full data sample collected with the Belle 
detector at the KERB asymmetric-energy e+e” (3.5 on 
8.0 GeV) collider [24], operating at the T(4S') resonance. 
The data sample consists of an integrated luminosity of 
711 fb“^, which corresponds to (771.6 ± 10.6) x lO^BB 
pairs. 

The Belle detector is a large-solid-angle magnetic spec¬ 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero¬ 
gel threshold Cherenkov counters (ACC), a barrel-like 
arrangement of time-of-flight (TOE) scintillation coun¬ 
ters, and an electromagnetic calorimeter (ECL) com¬ 
prised of CsI(Tl) crystals located inside a superconduct¬ 
ing solenoid coil that provides a 1.5 T magnetic held. 
An iron hux-return yoke located outside of the coil is in¬ 
strumented to detect AT® mesons and to identify muons 
(KLM). Two inner detector conhgurations were used. A 
2.0 cm beampipe and a 3-layer SVD were used for the hrst 
sample of 152 x IQ^BB pairs, while a 1.5 cm beampipe, a 
4-layer silicon detector and a small-cell inner drift cham¬ 
ber were used to record the remaining 620 x IQ^BB 
pairs |2S| . The detector is described in detail in Ref. |50| . 


The study was performed as a blind analysis based on 
simulated data. Monte Carlo (MC) samples were gener¬ 
ated with EvtCen m and the detector simulation was 
performed by GEANT3 [55]. Recorded beam background 
was added to the MC samples. The expected non-beam 
background is estimated using MC samples that describe 
all physics processes at Belle. A resonant T(4S') event at 
Belle produces a BB pair. Two samples of 6 —)■ c decays 
for B^B^ and B'^B~ events, respectively, each contain 
10 times the integrated luminosity of the data sample. 
Semileptonic b ^ u decays are simulated in a sample con¬ 
taining 20 times the integrated luminosity. Rare 5 —>■ s 
and other rare decays are described in another sample 
corresponding to 50 times the integrated luminosity of 
the data. Continuum e+e” qq (q = u, d, s, c) was gen¬ 
erated with PYTHIA [521 and included in the analysis in 
an MC sample containing hve times the integrated lumi¬ 
nosity of the data sample. Additionally, a high statistics 
sample of —?> X^tv containing 24 x 10® events was 

generated with a phase-space and ISGW2 [30] model. 

The signal MC sample is generated using BCL results 
for the vector and scalar form factors [5]. A total of 
84 X 10® R®R® events were generated with one meson 
decaying into the signal final state and the other decaying 
generically. 

No constraints on the r decay were applied. The signal 
MC sample corresponds to approximately 2000 times the 
expected B (R® —>■ Ti~T'^Vr) = 9.35 x 10“®. 

III. EVENT SELECTION 

The complete reconstruction of the B meson decay into 
the signal hnal state (Rsig) is not possible due to the pres¬ 
ence of at least two neutrinos. However, since the initial 
state of the e+e“ collision is completely defined by the 
momenta of the colliding leptons, we can constrain the 
signal side by fully reconstructing the other B meson 
(-Stag) in hadronic decay modes. Tracks and clusters in 
the event that are not assigned to the Stag after the suc¬ 
cessful reconstruction are assumed to originate from Rsig- 


A. Tag side 

This analysis uses the Belle hadronic full- 
reconstruction algorithm |31] based on the artificial 
Neural Network package NeuroBayes [32]. Neural net¬ 
works were trained to reconstruct R® and B^ candidates 
from a total of 1104 decay channels. 

Additional event shape variables are added to sup¬ 
press continuum events. R mesons of resonant events are 
nearly at rest in the center-of-mass (c.m.) frame, leading 
to a spherical distribution of their decay products. Con¬ 
tinuum events, on the other hand, produce back-to-back 
jets due the large available kinetic energy. Useful ob¬ 
servables that differentiate between the two event types 
are the thrust axis of the Rtag meson [33] and modified 
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Fox-Wolfram moments [31]. For this analysis, the thrust 
axis and the second modified Fox-Wolfram moment are 
included in the neural network for the full hadronic re¬ 
construction. If the algorithm does not succeed in recon¬ 
structing a 5° candidate, the event is discarded. 

Differences in the full reconstruction efficiency between 
MC and data, depending on the network output and Stag 
reconstruction channel, are observed [121 but depend on 
the tag-side reconstruction only; a correction factor is de¬ 
termined from charmed semileptonic decays of the signal- 
side B meson. 

The beam-energy-constrained mass, 

Mbc = \l /c^ 

is required to be greater than 5.27 GeV/c^, where Sbeam 
and pstag denote the beam energy and reconstructed 
three-momentum, respectively, of the Stag, evaluated in 
the e“'"e“ c.m. frame. With this requirement and the 
correction factor applied, we estimate a reconstruction 
efficiency of 0.18% from the signal MC sample, which 
is in very good agreement with the reconstruction effi¬ 
ciency of S° mesons in the Belle data sample [3T|. The 
neural network output, Ot|g S [0,1], is a continuous vari¬ 
able whose high (low) values correspond to candidates 
which are likely (unlikely) to be a true S meson. It is 
used at a later selection stage, as described below. The 
distributions of Mbc and ln(o5r®% for the three reconstruc¬ 
tion channels are shown in Figd^ where the green (solid) 
contribution shows the dominant charmed B —>■ X^. back¬ 
ground decays. No further requirements are applied to 
the Mbc distributions, while Mbc > 5.27 GeV/c^ is re¬ 
quired for the ln(ot®g) distributions. 

B. Signal side 

Only one-prong decays of the r lepton are considered 
in this search. For of a correctly reconstructed i?tag, 
there should be exactly two remaining oppositely-charged 
tracks in the detector. Additionally, the event should 
contain undetected (missing) momentum. Since the ini¬ 
tial state of the e+e” collisions is given by the four- 
momenta of the colliding leptons, the undetected mo¬ 
mentum can be measured. The missing momentum is 
defined as 

Pniiss — 2pbeam PStag PB^is ? 

where 2pbeam = Pe+ + Pe- is twice the beam momen¬ 
tum and Byis denotes the visible part of the Bgig meson. 
Tracks with low transverse momentum pt can curl in the 
solenoidal field and be detected as two tracks with oppo¬ 
site charge. Any two tracks with pt < 275 MeV/c with an 
angle between the momentum vectors, calculated at the 
points closest to the interaction point, below 15° and to¬ 
tal momentum difference less than 100 MeV/c are there¬ 
fore counted as a single track. We reduce the number of 


poor quality tracks by requiring that |dr| < 2.0 cm and 
\dz\ < 4.0 cm, where \dz\ and |dr| are the distances of 
closest approach of a track to the interaction point along 
the 2 :-axis and in the transverse plane, respectively. 

Electron identification [35| is performed by calculating 
a likelihood ratio using the matching of charged tracks 
with the shower position in the EGL, the shower shape, 
the ratio of the energy deposited in the EGL and the mea¬ 
sured momentum, the energy loss dE/dx in the CDG, 
and the Cherenkov light production in the ACC. Muon 
identification [32] is also done by evaluating a likelihood 
ratio. Clusters in the KLM are matched to charged tracks 
by extrapolation. Eor matched tracks, the difference be¬ 
tween expected and measured penetration depth and the 
transverse deviation of all KLM hits associated with the 
track are used in this likelihood ratio. Eor a charged 
track not identified as a lepton, a kaon veto is applied 
using a likelihood ratio that discriminates between kaons 
and pions m- The ratio is formed from the energy 
loss dE/dx in the GDC, flight time information from the 
TOF, and photon yield in the ACC. All remaining tracks 
are identified as a pion. Neutral pions are reconstructed 
from pairs of photons. The absolute difference between 
the invariant mass of the candidate and the nominal 
7r° mass, normalized to its uncertainty, must be below 
3.0. Photons are required to have energies in the lab¬ 
oratory frame greater than 50 MeV for the ECL barrel 
and 100(150) MeV for the forward (backward) endcap. 
Neutral pion candidates with at least one photon being 
used in the tag-side full reconstruction are discarded. 

Events are required to have exactly two oppositely 
charged particles within the allowed impact parameter 
range, with one additional track allowed outside the 
range. At least one charged pion is required. If the 
event contains two charged pions and neutral pion can¬ 
didates, we search for candidates. The charged pion 
with the lower momentum in the c.m. frame is combined 
with each neutral pion candidate and a mass vertex fit 
is performed. A pair that can be successfully fitted with 
< 20 is accepted as a charged meson if its in¬ 
variant mass is between 625 and 925 MeV/c^. If mul¬ 
tiple candidates are found, the candidate with a mass 
closest to the nominal p+ mass [53] is selected. Due to 
the broad mass range, not all t~ —>■ p~Vr events are 
correctly reconstructed. These events contain two oppo¬ 
sitely charged pions and a neutral pion in the final state 
and are miscategorized in the t~ —> ■k~Vt channel. Each 
event is reconstructed in one of the four reconstruction 
channels. In many t~ —>■ events, the momentum 

of the muon is too low to reach the KLM and thus is not 
identified as a muon. In most of these cases, the muon 
is identified as a pion so that the event is placed in the 
T~ — > 'k~Vt sample. 

Since Kl mesons do not completely deposit their en¬ 
ergy in the detector, charmed B decays with subsequent 
decays D or D —>• have the signal’s 

missing-momentum signature. A candidate is iden¬ 
tified as a cluster in the KLM without an associated 
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FIG. 1: Distributions of ln(o{r|g)(top) and Mbc (bottom) of the tag side candidate for the three reconstruction 
channels r —>• e (left), r —>■ tt (middle) and r —>• p (right) for signal and background, produced from Belle simulation. 
A requirement of Mbc > 5.27 GeV/c^ is applied in the left plot but no further requirements are applied other than 
successful reconstruction into one of the three modes. The histograms are produced from MG samples, normalized 
to the data sample size. The signal histogram is scaled by an arbitrary factor of 2000 for better visibility. The 
dashed vertical line indicates the minimal value required in the final selection; see Section 


HIE 


charged track. An ECL cluster without an associated 
charged track is associated with the Kl cluster in the 
KLM if it lies along the flight path extrapolated from 
the interaction point to the KLM cluster. As described 
below, the extra energy in the EGL is used to determine 
the signal yield. Therefore, only events with a Kj^ with¬ 
out energy deposition in the ECL are vetoed. 


C. Extra energy 


We extract the signal yield from a fit to the distribution 
of the energy deposited in the ECL (Eecl) by particles 
not used in the full reconstruction or by the two remain¬ 
ing charged signal tracks. To reduce background, the 
aforementioned photon energy requirements are applied. 
For signal decays, there is no additional energy deposi¬ 
tion, so Eecl peaks strongly at zero. Misreconstructions 
of i?tag lead to a small tail towards higher energy de¬ 
positions for true signal events. In contrast, most back¬ 
ground decays exhibit non-vanishing extra energy due to 
the presence of additional neutral particles. 


D. Boosted decision trees 

Final event selection uses requirements on three vari¬ 
ables: ln(ojr®g); missing mass squared (M^j^^ = p^igg); 
and the output of the boosted decision tree (BDT). For 
each T reconstruction mode, one BDT is trained using the 
TMVA framework [55]. All use different input variables, 
background training samples, and BDT growth parame¬ 
ters. The signal training sample consists of 3 x 10^ events 
out of the complete signal MG sample. To improve the 
training, events are required to have Eecl < 1 GeV and 
Stag is required to have a quality of ln(o5;®g) > —7. One 
additional track outside the impact parameter require¬ 
ment is allowed. 

Another 3 x 10^ signal events are used for performance 
tests of the BDT for receiver-operation characteristics 
(ROC) calculation and overtraining evaluation. 

The input variables of the BDT used in the t~ — >■ 
e~VfyT selection are the magnitude of the three-momenta 
of the pion and electron, the squared lepton-pair momen¬ 
tum transfer and different combinations of all 

available four-momenta. The momentum transfer can be 
calculated using the fact that both B mesons are at rest 
in the c.m. frame, which implies = —PStag 

q = ps^ig — p,r. Due to the low efficiency of the full re¬ 
construction, an additional signal sample is used with 
2 X 107 r- ^ e- Vey-r events on the signal side. The 
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background training sample consists of charmed de¬ 
cays and —)■ Xuii'e decays. The input variables are 

linearly decorrelated before their use in the training. 

The background sample used in the t~ — >■ se¬ 

lection BDT contains b ^ c decays and semileptonic 
b ^ u decays of B^ mesons. Principal-component analy¬ 
sis (PCA) [38] is applied to the input variables. PCA is 
an orthogonal transformation which rotates a sample of 
data points such that the variability along the new axis is 
maximized. In this way, the variables are decorrelated. 
The input variables are the missing energy, < 7 ^, 

the absolute three-momentum of Bvis in the c.m. frame, 
combinations of available four-momenta, and the number 
of unused neutral pions in the event. 

The BDT training for the t~ —>■ p~ Vr selection uses the 
same sample size of 6 —>■ c decays, but not semileptonic 
b ^ u decays. The correlation of the input variables is 
again reduced by a PCA transformation. The variables 
used in the training are the missing energy, 

and combinations of the available four-momenta in the 
decay. 

The performance of the three final BDTs is shown in 
form of the ROC curves in Fig. 


Background rejection versus Signal efficiency 



FIG. 2: Background rejection versus signal efficiency 
determined on the testing sample, from Belle 
simulation. 


E. Final selection 

The final selection criteria are determined from MC 
samples by maximizing individually the expected signifi¬ 
cance of each single r reconstruction mode. We perform 
a scan over three variables simultaneously to obtain the 
optimal selection: ln(ot|g), M^^gg, and the BDT output. 
We require ln(ot|g) > —7 for the leptonic t reconstruc¬ 
tion and ln(o5r|g) > —5 for the hadronic r reconstruction, 
as shown in Fig.[^ A minimum requirement on the Af^jgg 
is applied to reject semileptonic B — > t:(,v events, which 
have the same final state as the signal decay: since no 


energy is deposited in the ECL, decays of this type peak 
at zero extra energy. Also, as there is only a single neu¬ 
trino in these decays, M^jgg peaks at zero, unlike the 
case for signal decays, which contain at least two neu¬ 
trinos. We require AfAgg to be greater than 2.2 GeV/c^ 
in the electron channel, 0.0 GeV/c"^ in the pion channel, 
and 0.6 GeV/c^ in the p channel. The BDT output is 
the last variable used in the scan. The expected signif¬ 
icance is calculated as a/— 21n(£o/>Ci); the likelihood is 
given by Ck = Hbins Poi®son(nobs|»T'bg + k ■ risig), where 
f^obs = ■'T'bg + nsig is the best estimate from the MC sam¬ 
ples. 

The efficiency of the final selection is determined from 
MC to be 4.57 x 10“'^. The dominant reconstructed r 
decay modes and their relative occurrences are listed in 
Table H 

TABLE I: Signal reconstruction by r decay modes. 

Percentages are obtained from signal MC and sum to 

100 %. 


T decay 

Relative Occurrence (%) 

P~Vt 

29.54 

e~Ve^T 

29.43 

'K~Pt 

16.70 


13.21 

a^Vr 

8.72 

other 

2.4 


The dominant background in the Iow-Aecl region 
arises from 5° —and B^ —>■ D^*^p decays with a 
subsequent decay of D —>• K^tt. The is undetected in 
these cases and the resulting decay signature resembles 
that of the signal. No explicit selection is available to 
further suppress decays of this type. 


IV. SYSTEMATIC UNCERTAINTIES 

In the computation of the significance level and upper 
limit, systematic uncertainties are included in the like¬ 
lihood as nuisance parameters. The likelihood is built 
from probability density functions (PDFs) determined 
from MC predictions of each background sample, as de¬ 
scribed in Section [V] All systematic uncertainties are as¬ 
sumed to be Gaussian-distributed and are evaluated at 
one standard deviation (u). 

Uncertainties of the particle identification and of the 
correction factor needed for the full reconstruction effi¬ 
ciency are included as a flat effect over all bins in Aecl- 
All other uncertainties are included in a bin-by-bin fash¬ 
ion. A constant uncertainty of 0.35% has been deter¬ 
mined for each charged track with pt > 0.2 GeV/c. 
Tracks below that threshold have to be treated differently 
depending on the track momentum |39j . The uncertainty 
on the number of produced R-meson pairs is 1.4%. The 
uncertainty due to the veto is determined by vary- 
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ing the efficiency by its uncertainty. The branch¬ 
ing fractions of the dominant backgrounds are varied by 
their errors stated in Ref. [53] to determine the effect 
on the MC prediction. The uncertainty on the correc¬ 
tion factor of the tagside reconstruction is determined in 
Ref. [19] and applied to the samples. The discrepancy 
between inclusive and exclusive \Vub\ measurements has 
been included as a flat but asymmetric uncertainty in 
the B —>■ Xu(.v sample of {1li^)%. An uncertainty of 
±10% is applied to the branching fractions in the MC 
sample of rare b ^ s and other rare B decays. Addi¬ 
tionally, decays of type B —>■ X^tv are present in the 
final event selection. The contribution to the Aecl dis¬ 
tribution is evaluated from the MC sample assuming a 
B{B^ — )■ = 1.5 X 10“^ and found to be small; a 

relative uncertainty of ±50% is applied. Statistical un¬ 
certainties in the PDF shape due to finite MC sample 
size are included in a way similar to the approach by 
Barlow and Beeston [30]. Instead of using one Poisson 
constraint per background sample per bin per r decay 
channel, only one constraint term per bin per channel is 
used. The uncertainty introduced by this approximation 
is negligible for bins with non-vanishing content and re¬ 
duces the amount of computation time needed. Instead 
of the finite MC uncertainty, the fit error is included as 
a systematic uncertainty for the dominant 6 —>■ c contri¬ 
bution. The theoretical uncertainties of the signal form 
factors /+ and /° are included by generating additional 
signal MC with one form factor fixed and the other varied 
by its Icr uncertainty. The relative uncertainties deter¬ 
mined in this way are combined into a single uncertainty 
estimate. The systematic uncertainties due to the track¬ 
ing efficiency and particle identification affect only the 
overall efficiency and are only included in the calculation 
of the upper limit. The relative effect on the branch¬ 
ing fraction is determined by repeatedly fitting modified 
PDFs to data. The PDFs are modified by replacing each 
background contribution with the respective contribution 
where the systematic effect is applied. For each system¬ 
atic uncertainty, two fits are performed for the positive 
and negative deviation. The maximum, absolute devia¬ 
tion is quoted in Table |n| 


V. RESULT 

A binned maximum likelihood fit is performed to Aecl 
in bins of 0.15 GeV. Due to similar shapes in the back¬ 
ground predictions, all background contributions except 
for the dominant 5 —>■ c transitions are fixed to the MC 
prediction. Possible errors introduced by this approach 
are accounted for as systematic uncertainties. The fit 
is performed simultaneously in all three reconstruction 
modes. The signal strength parameter p, is constrained 
between the three modes while the background contribu¬ 
tions of the three reconstruction modes are floating pa¬ 
rameters. The fit result of the B^ —)■ X^ background con¬ 
tribution agrees well with the prediction obtained from 


TABLE II: Effects of the systematic uncertainties on 
the branching fraction. 


Source 

Relative error (%) 

Particle ID 

2.4 

Track efficiency 

0.7 

N{BB) 

1.4 

Kl veto 

3.2 

BG B 

2.8 

Tagside 

4.6 

\Vub\ 

2.8 

Rare processes 

2.0 

B — >■ XuTU 

2.2 

Background fit 

0.2 

Signal model 

1.8 

Total 

8.3 

the MC sample. The signal strength has been chosen 
such that p = 1.0 corresponds to B{B^ —Tr^r+i/,-) = 
1.0 X 10“"^. We obtain a best fit of ^ = 1.52±0.72, corre¬ 
sponding to 51.9 ± 24.3 signal events. The fit results by 

T reconstruction mode are 

listed in Table |lll| The -Eecl 

TABLE III: Fit results for signal yield. Total and split 

by r reconstruction mode. 

Mode 

Signal Yield 

e 

13.2 ±6.2 

TT 

30.6 ± 14.3 

p 

8.1 ±3.8 

Total 

51.9 ±24.3 


distribution and fit results are shown in Fig.[^ 

The significance of the measurement is obtained from 
a pseudo MC study. A test statistic based on the profile 
likelihood ratio is used. The likelihood is built in bins of 
0.15 GeV in Eecl- The binned likelihood is given by 

C = n n Poisson(nc6|r'cb) • fp{.ap\ap), (3) 

c b pGP 

where the indices c and b label the reconstruction chan¬ 
nel and bin in Eecl, respectively, and P denotes the set 
of systematic uncertainties p that are included as nui¬ 
sance parameters ap in the calculation of the number of 
expected events Vcb per channel per bin. The nuisance 
parameters are parametrized as a relative effect on the 
nominal template prediction, assumed to be Gaussian- 
distributed with the nominal value being the global ob¬ 
servable ap. The number of events in the background- 
only hypothesis is determined from MC simulation and 
a fit to data for the dominant b ^ c background. The 
likelihood is constructed using the HistFactory tool in 
the RooStats package [4T| |42| . 

The distribution of the test statistic is obtained by 
pseudo-experiments. A full frequentist approach is used 
in both the computations of the significance level and 
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FIG. 3: Distributions of F'ecl in the three r reconstruction modes. The signal and b ^ c contributions are scaled 

according to the fit result. 


the upper limit. First, the likelihood is fitted to data to 
obtain the maximum likelihood estimates (MLEs) of all 
nuisance parameters on data. In each pseudo-experiment 
generation, the nuisance parameters are fixed to their 
respective MLE. In the subsequent maximization of the 
likelihood, the nuisance parameters are free parameters. 
The global observables are randomized in each pseudo¬ 
experiment. 

Using pseudo-experiments, the p-value of the 
background-only hypothesis for data is determined 
and the significance level Z is computed in terms of 
standard deviations as 

where is the cumulative distribution function of the 
standard normal Gaussian. 

We observe a signal significance of 2.8a, not includ¬ 
ing systematic uncertainties in the calculation. Including 
all relevant systematic effects results in a significance of 
2.4cr. For this result, the test statistic has been computed 
on 10 000 background-only pseudo-experiments. 


Given the level of significance of these results, we invert 
the hypothesis test and compute an upper limit on the 
branching fraction, pseudo-experiments are generated 
for different signal strength parameters for both signal- 
plus-background and background-only hypotheses in or¬ 
der to obtain CLg+b and CLf,, respectively. The upper 
limit is then computed using CLs = CLg+b/CLi, [iH] . 
where a scan over reasonable signal strength parame¬ 
ter values is performed. At each step, 10 000 pseudo¬ 
experiments have been evaluated for both hypotheses. 


At the 90% confidence level, we obtain an upper 
limit of B (i?° —>■ < 2.5 x 10“^. The upper 

limit at the 95% confidence level has been computed to 
B (i?° —)■ < 2.8 X 10“^. This result is the first 

result on B {B^ —>■ and is in good agreement 

with the SM prediction. 
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